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一 种 适用 于 烷烃 高 压 小 相 藉 度 的 推算 模型 
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摘 要 本 文 基于 绝对 速率 理论 和 SRK 状态 方程 , 建立 了 一 种 可 以 精确 计算 烧 烃 高 压 液 相 黏度 的 推算 模型 .选取 10 种 
正 构 烷 烃 (Ci-Cio) 的 液 相 黏 度 作为 研究 对 象 ， 温度 范围 为 90~448 K， 压 力 范 围 为 0.01~254.4 MPa, 利用 本 文 所 提出 的 
黏度 模型 对 其 进行 了 计算 . 计算 结果 表明 : 本 文 的 黏度 推算 模型 对 这 10 种 烷烃 禁 度 的 计算 值 与 实验 值 之 间 的 平均 绝对 偏 
差 小 于 3%, 验证 了 模型 的 精确 性 . 
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A. Viscosity Model of Liquid n-Alkanes at High Pressures 


LIU Xiang-Yang ZHU Chen-Yang Mao-Gang 
(MOE Key Laboratory of Thermal-Fluid Science and Engineeri ian Jiaotong University, Xi'an 710049, China) 
Abstract In this work, a new viscosity mode eh had a high accuracy of liquid n-alkanes 
at high pressures was presented. The we with the absolute rate theory and the SRK 


equation was applied to correlate the ps 
90-440 K at pressures up to 254.4 d The overall average absolute deviation between exper- 
imental data and calculated values ss than 396, which verified the precision of the present model. 


les of ten n-alkanes in the temperature range of 
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LEM d 
H Y 
黏度 是 流体 基本 的 迁移 性 质 之 一 ， 对 工业 生产 通过 构建 其 T —(— P 
及 过 程 优化 方面 都 有 着 最 重要 的 意义 . 烷烃 作为 一 THIERA CAE EREE 
种 常见 的 民用 燃料 与 化 工 原料 ， 其 液 相 黏度 数据 在 J T jRSEHESRUR SUE EHE 
油田 开采 、 化 工 设计 等 领域 都 有 着 重要 的 应 用 . HO M Men 
Wr, ERRER EERIE EIE T Vo c (D ffartins FRI n BETERLPUS 


量 工作 "但 是 实验 研究 只 能 提供 离散 点 , WAS ^—— 对 于 牛顿 流体 ,根据 Eyring 的 绝对 速率 理论 ， 


真正 满足 实际 需求 基于 实验 数据 ， 开 发 高 ; 可 得 其 务 度 的 表达 式 为 
理论 计算 方法 是 获取 务 度 数据 的 有 效 途 径 . 

虽然 现在 国际 上 学 者 们 提出 了 许多 液态 流体 区 )= 生 加 ( 徐 ) a) 
度 的 理论 计算 方法 ,但 大 多 只 适用 于 常 压 流体 . 对 d RT 


Tum HRULURBUSRHE, BUNIUJBQERGETRSUES 2000 年 
Quiiones-Cisneros 等 [1] dH f$ EE SER C s EA 
和 2001 4E Allal 等 !? 提出 的 自由 体积 理论 黏度 模 
型 ， 此 后 ，Martins 等 03 基于 绝对 速率 理论 和 状 


AF, 7 为 动力 黏度 , h 为 普 朗 克 常 数 ,N AARI 
德 罗 常 数 , V 为 流体 的 摩尔 体积 ，AF* 为 流动 过 程 
中 的 摩尔 活化 能 , RR 为 气体 常数 , T 为 热力 学 温度 . 


态 方程 ， 建 立 了 一 种 可 适用 于 高 压 液体 黏度 的 理论 人 
推算 模型 ， 对 正 构 烷 烃 务 度 计 算 的 平均 绝对 偏差 为 hN AF; 
1.22%。 然 而 ， 对 于 每 个 温度 ， 需 要 通过 实验 数据 拟 m= re (FE) Q2) 


合 出 一 套 模 型 参数 , 大 大 限制 了 该 模型 的 实际 应 用 . 
本 文 以 Martins 等 013] 的 务 度 理论 模型 为 基础 ， 其 中 , 下 标 0 表示 参考 状态 . 
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为 消去 式 (1) 中 的 常数 , 将 式 (1) 与 式 (2) 相 比 
并 整理 可 得 
n= m exp (1) (3) 


假设 流体 的 摩尔 活化 能 之 差 AF" — AF; 等 于 其 摩 
尔 自由 能 之 差 , 即 
AF*—-AF;—-AF-F-H (4) 


将 式 (4) 代入 式 (3) 中 , 并 取 该 模型 下 流体 的 自由 能 
为 Helmholtz 自由 能 A, 可 得 
Vo AA 

n = Tho Exp (27) (5) 
AB, AA 为 流体 给 定 状态 与 参考 状态 下 的 Helmholtz 
自由 能 之 差 . 对 式 (5) 中 的 摩尔 体积 比 V/V 做 维 里 
形式 的 展开 , 则 有 

Vo 


T = 1+ B)pt Bi)? (6) 


AF, 2 为 压力 , B(T) 和 B(T) 为 只 与 流体 温度 有 


Chung 提出 的 关联 式 得 到 (141 
m0 = 4078.57* (mm " (7) 
AF, mo HERABE, KISTE, T 为 热力 
学 温度 ， 友 为 临界 摩尔 体积 BY 为 碰撞 积分 ,为 
与 分 子 形状 和 极 性 有 关 的 因子 , 其 表达 式 为 
F, = 1 — 0.2756w + 0.05903544 + K (8) 


其 中 , w 为 偏心 因子 ，# 为 对 高 极 性 分 子 的 修正 项 ， 
对 于 非 极 性 流体 可 取 0, p HAERE, 可 由 下 


式 给 出 
H 
Lr = 131.3 — = x 
(ome X 


式 中 , n WARE. 碰撞 积 分 O, 取 Lennard-Jones 
势能 函数 , 则 可 由 下 式 计 算得 到 US 
R = A/T*P + C/exp (DT*) + E/ 


exp (FT*) + RT* (ST*W — P) un 


其 中 , 4 = 1.16145, B = 0.14874, C = 0.52487, D = 
0.727320, E = 2.16178, F = 243787, R = —6435 x 
10-4, S = 18.0323, W = —0.76830, P = 7.27371, 
T* = 1.2593T,, 
将 式 (6), (7) 代入 式 (5) 中 得 
(MT)? 


Fe 
f= .065 e 


R 
BT) BC] x exp (25) 


à | AY 
选取 参考 状态 为 理想 状态 , 则 式 (5) 中 gie 


AF, AP 为 剩余 摩尔 自由 能 ,可 通过 选取 一 定 的 状 
态 方 程 求 得 . 

对 于 立方 型 的 状态 方程 ， 可 写 为 如 下 的 通用 
形式 

RT a 
P—V-b V2} ubV wb uo) 

将 压缩 因子 Z= PV/RT 代入 上 式 中 ， 可 将 状态 方 
程 改写 为 如 下 形式 


Z? — (1+ B* — uB*) Z? + (A* + wB?— 


wB*? 


uB* -uB?) Z — A*B* wB* =0 (13) 


者 选取 状态 方程 为 SRK 方程 , ERR u= 1, w= 0, 
则 式 (13) AEA 


eem + (A*— B* — B®?) Z — A*B*=0 (14) 


A* = ap/ PT? (15) 
B* =bp/RT (16) 


a = 0.42724 R? T2 1 +f (1 = | /fr (17) 


b = 0.08664 RT./ pe (18) 
cS 
其 中 X_。 
N fo = 0.48 + 1.574w — 0.1764? (19) 


N 状态 方程 为 SRK 方程 , 则 剩余 Helmholtz H 
能 为 


dc. 


将 式 (15). (16) 代入 式 20) 中 可 得 


ZRT 


AR B a/b RT RT 
RT ^ "A[ZRT bp 


" nis) e 


将 上 式 代入 到 式 (11) P, 则 可 得 实际 流体 黏度 的 推 


F. (MT)? 
n = 4.0785— 3 7.——— [1 + Bi(T)p + Ba(T)p?] x 
Ve Qy 
RT ZRP YOPE (22) 
ZRT — bp \ ZRT + bp 
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表 1 本 文 使 用 的 黏度 数据 
Table 1 The used viscosity data 


Liquid T/K p/MPa aNp PRef. 
CH4 100-140 0.1~50 43 1 
C2H6 100~210 0.1~60 126 2 
C3H8 90~240 0.01~100 281 3 
n-CaHio 280~348 2-70 105 4], [5] 
n-CsHi2 303.15~380 0.1~251.6 118 4]~[7] 
mn-C6Hi4 303.15-.420 0.177250 125 [4], [5], [7] 
n-C7H16 300~360 0.1~100 84 5], [8] 
n-CsHig 303.15~448 0.1~253.1 149 [4], [5], [7] 
n-C9H20 300~420 0.1~50 221 
n-C10H22 280~440 0.1~254.4 312 [5]. [7], [9], 10] 


^Np— 拟 合 数据 点 数量 ; "Ref. 一 数据 文献 来 源 . 


表 2 Martins 等 模型 参数 
Table 2 The parameters in VY etis model 


ZA 


Liquid T/K Bi B3 DEREN T/K Bı B2 
CH4 100 4.60E-04 5.81E-12 WF 7Hi6 340 4.19E-04 4.17TE-12 
110 1.26E-04 1.67E-12 343 3.69E-04 2.49E-12 
120 4.23E-05 o aR 360 1.67E-04 1.75E-12 
130 1.69E-05 D n-CsHis 303 1.34E-02 1.18E-10 
140 7.7AE-06 E3 320 4.41E-03 4.29E-11 
C2H6 100 1.02E+01 y. PÕE-07 323 3.76E-03 3.60E-11 
110 9.93E-01 Sd 8.78E_09 340 1.44E-03 1.43E-11 
m. 120 LAGE] T" 1.13E-09 348 9.A7E-04 9.29E-12 
qu 130 2. 2.09E-10 360 5.53E-04 5.92E-12 
= 140 | 5.08E-11 373 3.03E-04 2.85E-12 
150 w 0E-03 1.54E-11 3 2.18E-04 2.58E-12 
160 8.25E-04 5.5TE-12 GN 1.18E-04 8.91E-13 
170 3.37E-04 2.29E-12 * 423 5.29E-05 3.08E-13 
CoHe 180 1.52E-04 1.07E-12 | n-CsHis\ < 448 2.75E-05 8.89E-14 
190 7.43E-05 5.53E-13 "eg 300 7.48E-02 7.72E-10 
C 200 4.15E-05 3.19E-13 . es 310 3.51E-02 3.50E-10 
210 2.20E-05 1.83E-13 320 1.75E-02 1.67E-10 
: CsHg 90 1.03E--06 149E-02 . 330 9.14E-03 8.52E-11 
á 100 1.89E+04 N 340 4.99E-03 4.61E-11 
110 7.80E--02 PRO 350 2.84E-03 2.62E-11 
120 5.72E--01 6.30K-07 360 1.67E-03 1.55E-11 
130 6.49E--00 6.46E-08 370 1.02E-03 9.37E-12 
140 1.03E4-00 9.43E-09 380 6.36E-04 5.84E-12 
150 2.14E-01 1.81E-09 390 4.12E-04 3.65E-12 
200 9.88E-04 7.03E-12 400 2.73E-04 2.38E-12 
240 7.1TE-05 5.10E-13 410 1.85E-04 1.61E-12 
n-C4H10 280 7.89E-05 7.76E-13 420 1.28E-04 1.13E-12 
300 3.27E-05 2.54E-13 n-C19Ha3 280 2.08E-4-00 2.38E-08 
340 7.45E-06 7.39E-14 290 7.89E-01 8.75E-09 
320 1.50E-05 1.29E-13 297 3.84E-01 3.71E-09 
323 1.53E-05 8.54E-14 300 3.22E-01 3.37E-09 
348 6.40E-06 5.08E-14 303 2.20E-01 2.85E-09 
n-C5Hi2 303 1.54E-04 8.39E-13 310 1.40E-01 1.43E-09 
313 9.39E-05 6.33E-13 311 1.28E-01 1.52E-09 
318 7.45E-05 5.91E-13 313 1.12E-01 1.10E-09 
320 6.76E-05 5.96E-13 320 6.48E-02 6.75E-10 
323 6.17E-05 4.72E-13 323 4.87E-02 5.17E-10 
328 5.07E-05 3.99E-13 328 3.72E-02 3.80E-10 
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续 表 
Liquid T/K Bi Bə | Liquid T/K Bı Bə 
340 3.10E-05 2.77E-13 330 3.18E-02 3.30E-10 
360 1.52E-05 1.54E-13 340 1.65E-02 1.68E-10 
380 7.99E-06 9.37E-14 343 1.36E-02 1.41E_10 
mn-C6Hi4 303 7.22E-04 4.39E-12 348 1.00E-02 7.65E-11 
313 4.06E-04 3.42E-12 350 8.92E-03 8.95E-11 
323 2.53E-04 2.03E-12 358 5.41E-03 6.16E-11 
348 8.35E-05 7.81E-13 360 5.00E-03 5.09E-11 
373 3.44E-05 2.88E-13 370 2.90E-03 3.01E-11 
380 2.58E-05 3.51E-13 373 2.50E-03 1.96E-11 
400 1.36E-05 1.91E-13 380 1.74E-03 1.83E-11 
420 7.59E-06 1135-13 390 1.07E-03 1.15E-11 
n-CzHis 300 3.69E-03 3.77E-11 400 6.82E-04 7.45E-12 
303 3.09E-03 2.43E-11 420 2.97E-04 3.32E-12 
320 1.15E-03 137E-11 I A40 1.40E-04 1.65E-12 
323 9.91E-04 7.14E-12 c 


€— má 2 为 待定 参数 . 将 式 (23) 代 
AR (22) 中 , 可 得 
利用 Martins 等 19) 提出 的 务 度 模型 对 D 


构 烷 烃 (C1 ~C10) FERAE F 8928 BERE 联 F, (MT)? . B 
共计 1564 个 数据 点 , 温度 范围 为 9 人 poto T [1 + pexp (ea + AT; 7*) + 


Ul Un. 拟 合 
图 1~2 所 示 。 


围 为 0.01~254.4 MPa, RJ Hf? 
p? exp (o -- BT; ??)] " 


RR ZRT | 
um e 
ZR 4 bp VZRT b 
RR 1 F ;) 
xt (29 中 只 有 6 个 与 温度 和 压力 无 关 的 待定 参 


、ml=12), 相对 于 Martins 等 H3] 提出 的 
E 算 模型 ， 式 (24) 普 适 性 更 强 , 适用 于 不 同 温 


图 1 InB1 随 温 度 的 变化 曲线 


Fig. 1 Temperature dependence of lInB1 


通过 图 1~2 可 以 发 现 , 对 于 液 相 烷烃 ，Martins 

等 9 提出 的 黏度 推算 模型 参数 Bl 和 E 2 RUE 

的 增 大 呈现 递减 趋势 , 且 可 拟 合 为 对 比 温度 T; 的 关 5 
图 2 InB2 随 温度 的 变化 曲线 

nB =+ ?, i—1,2 (23) Fig. 2 Temperature dependence of InB» 


HI 
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表 3 本 文 模型 的 参数 值 及 计算 结果 


Table 3 The parameters and calculated results of our model 


n m aAAD/% 

.472 $ 
7.47 003 ET 
8.536 0.896 

a 1.1 
6.839 00 "m 
4.620 1.356 
6.572 1.168 " 
6.731 1.174 
10.345 0.864 "m 
0.694 4.386 

12i 1.2 

7.230 09 js 
—23.14 —0.514 

11.2 * 

5 0.897 m 

1.138 3.128 
11.516 0.928 1.87 
0. 3.325 


Liquid Qi 
三 1 —21. 
CH4 i 95 
{= —41.08 
= —21.02 
C2H6 ? 0 
j= —37.01 
CsHs T= —20.41 
C= —39.51 
n-C4Hio ve —24.29 
g= —32.21 
=] —21. 
n-C5 H12 : or 
i= —9.28 
= —25.1 
n-C6Hi4 i 5.13 
UE —31.76 
n-C7Hi6 i= —25.48 
i—2 一 30.77 
n-CsHis ys —17.36 
22 9.59 


"TT i= 1 —20.20 AÀ 


i—l1 i 
n-C10H22 qu 
i= 2 .91 


: 1.660 271 
44.54 —0.496 
1.064 1.335 0.96 
6.584 1.409 
.A11 1.32 
i 229 1.73 
5.044 1.635 


Np - 
^AAD = 5 [ricalj 一 velia /, zs 计算 值 与 实验 值 的 平均 绝对 偏差 


j= iS 


使 用 最 小 二 乘法 拟 合 得 到 了 适用 于 上 述 10 种 烷 
烃 的 待定 参数 oi. Bi. nli — 1,2), WK 3 所 示 , 在 此 
基础 上 ,利用 式 (24) 对 10 种 烷烃 的 液 相 黏度 进行 
了 计算 , 并 与 文献 中 的 实验 值 进行 了 比较 , 计算 结果 


A-74(77100 K) o rer(7=100 K) 
-e- n4 ( T7110 K) o ner(7=110 K) 
A-na( T7120 K) ^ mrer(T=120 K) 
nns T7130 K) v (77130 K) 
ns T7140 K)o ma(T-140 K) 


p/MPa 
图 3 CH4 黏度 计算 值 与 文献 值 比较 


Fig. 3 Comparison between the calculated values and the 


literature data for the viscosities of CH4 


如 表 3 和 图 3E. 结果 显示 , 式 (24) 对 10 种 
正 构 烷烃 液 邦 猎 度 的 计算 值 与 实验 值 吻合 良好 ， 两 
者 之 NU S SHUT 3%. 


é 


—— qa(T-100 K) 4 44(T7150 K) -®-Ncal(T=190 K) 
D qma(T-100K) 4 mQ4(T-150K) © Mret(T=190 K) 
—e- qa (T7110 K) >- noa(T=160 K) —9— nca(T=200 K) 
© qa(T-110 K) > ma4(T-160 K) + Mret(T=200 K) 


——1qm4(T-120 K) -9—5,(77170 K) =- Na(T=210 K) 
A MelT=120K) o mel(T=170K) x m4(T-210 K) 
Y- na( T7130 K) —— na(T=180 K) 
Y MelT=130 K) * ma4(T-180 K) 
—e— na(T=140 K) 


PEEL E 
ER 

"| 放电 
0 = m . 1 

0 10 20 30 40 50 60 


p/MPa 


图 4 C2He SEBETETEBLS SCIRE EGRE 
Fig. 4 Comparison between the calculated values and the 


literature data for the viscosities of C2 Hg 
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——54(T-90 K) —v—54(7T-120 K) ——74(77150 K) 
o ma(T-90K) Y pee(T=120 K) > mar(T=150 K) 
—e— Na(T=100 K) —9—5.,(T-130 K)-®-.a(T=200 K) 
o ma(T-100K) © m4(T-130 K) © Mret(T=200 K) 
A i(T7110 K) —7j,4(T7140 K)——74(1-240 K) 
^ mr(T=110 K) 4 nma4(T-140 K) * na ( 77240 K) 


o 


p/MPa 
图 5 CsHs 黏度 计算 值 与 文献 值 比 较 


Fig. 5 Comparison between the calculated values and the 
literature data for the viscosities of C3Hg 


A qa(1-280 K) ^. n4(T-280 K) 
—- qa(T-300 K) v. m 4(7-300 K) 
350 F —e- 4, (7-320 K) o. (7-320 K) 
—— qma(1-323 K)o n,4(7-323 K) 
300 - —— na (T7340 K) ad Nret(T=340 K) 
—e- qa (1-348 K)o nue(T=348 K) 


100 
p/MPa 
图 6 C4aHio A REIES OC EGRE 


Fig. 6 Comparison between the calculated values and the 
literature data for the viscosities of C4H1o 


450 
-A-na( T7313 K) 4,,(7-323 K) -© 7,4(7—360 K) 


400| ^ «(7-313 K) © ,,(7-323 K) © Mer(T=360 K) 
A8 54(77318 K) Y- NalT=328 K) #7.a(T=380 K) 
D qa(T-7318 K) Y (7-328 K) * Met(T=380 K) 
-Nea T7320 K) > Na(T=340 K) 
4 ia (T=320 K) > Me(T=340 K) 


350 


n/uPa-s 


0 10 20 30 40 50 60 70 80 
p/MPa 


图 7 CsHio 黏度 计算 值 与 文献 值 比较 


Fig. 7 Comparison between the calculated values and the 


literature data for the viscosities of C5 H12 
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图 8 CeHi4 黏度 计算 值 与 文献 值 比较 


Fig. 8 Comparison between the calculated values and the 
A data for the viscosities of CGH14 
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NSE 9 Comparison between the calculated values and the 
N literature data for the viscosities of C7H16 
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图 10 CsHis 黏度 计算 值 与 文献 值 比较 


Fig. 10. Comparison between the calculated values and the 


literature data for the viscosities of CsHig 
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1800 具有 较 高 计算 精度 。 
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图 11 CoH»o 黏度 计算 值 与 文献 值 比较 


Fig. 11 Comparison between the calculated values and the 
literature data for the viscosities of CogH20 
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图 12 CioHz» Ri BETETUELS CHE E EX 
Fig. 12 Comparison between the calculated values and the \ 
sA 


literature data for the viscosities of C1i0H22 
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3 结论 ~ 

本 文 将 绝对 速率 理论 和 SRK 状态 方程 相 结 合 ， 
并 通过 构建 模型 参数 与 温度 之 间 的 关联 函数 ， 建 立 
了 一 种 可 以 精确 计算 烷烃 高 压 液 相 黏度 的 推算 模型 . 
以 温度 为 90—448 K、 压 力 为 0.01~254.4 MPa 范围 
内 , 10 种 正 构 烷烃 (Ci Cio) 的 液 相 黏度 为 基础 数 
据 , 利用 最 小 二 乘法 获取 了 模型 参数 . S088 BETESY 
模型 对 这 10 种 烷烃 液 相 黏 度 的 计算 值 与 实验 值 之 间 
的 平均 绝对 偏差 小 于 376, 表明 本 文 的 黏度 推算 模型 
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